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"sj - . 1 INTRODUCTION 

o 

H ' Gamma-ray bursts (GRBs) are one of the brightest phenom- 
ena in the Universe. Although the majority of GRBs were 
detected at energies ranging from hundreds of keV to several 
MeV, they w ere also observed at much higher energies up to 
^ I tens of GeV (|Hurlev et al.lll994 ). The advent of the Fermi 
J-h . Large A rea Telescope (LAT) with its unprec edented sen- 
sitivity (|Atwood et al.ll2009l ; iBand et alj 120091 ') has greatly 
increased ou r capability to stud y the high-energy emission 
from GRBs l|Omodei et alJl2009l ). 

The high-energy (HE) emission was detected both 
in the prompt and afterglow phases of GRBs. However, 
its origin is still unclear: it could be produced in the 
internal/external shocks via leptonic or hadronic mech- 
anisms, or in the process of d i ssipat io n of the Poynt- 



ABSTRACT 

We search the Fermi-LAT photon database for an extended gamma-ray emission 
which could be associated with any of the 581 previously detected gamma-ray bursts 
(GRBs) visible to the Fermi-LAT. For this purpose we compare the number of photons 
with energies E > 100 MeV and E > 1 GeV which arrived in the first 1500 seconds 
after the burst from the same region, to the expected background. We require that the 
expected number of false detections does not exceed 0.05 for the entire search and find 
the high-energy emission in 19 bursts, four of which (GRB 081009, GRB 090720B, 
GRB 100911 and GRB 100728A) were previously unreported. The first three are 
detected at energies above 100 MeV, while the last one shows a statistically significant 
signal only above 1 GeV. 
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emission could shed light on the onset of GRB and its im- 
mediate interaction with surrounding interstellar medium. 

GRB observations in the HE band are one of the key 
science topics of Fermi-LAT. The Fermi Gamma-Ray Burst 
Monitor (GBM) can initiate autonomous slew of the space- 
craft to provide the best conditions for a dedicated GRB 
observation. Also, knowledge of time and position of a 
GRB (either provided by GBM or by other observations) 
makes it feasible to search for GRBs in the LAT photon 
database where they would manifest themselves as spatially 
and temporally compact clusters of photons (B and et al.l 
l2009h . There has been 20 LAT detections of GRBs as of 
the time of writing (Feb 2011). 

Observations of several very bright bursts (GRB 
080916C, GRB 090510, GRB 090902B, GRB 090926A) con- 
taining more than 100 photons with energies in excess of 
100 MeV allowed one to study the spectral prop erties of the 
HE em is sion and even their temporal ev olution (|Abdo et al.l 
l2009d lal: lAckermann et all |2010| . l201ll ) . In all these cases 
the HE emission demonstrated a delay (of several seconds 
for long bursts, tenth of a second for the short burst GRB 
090510) with respect to the prompt emission in the sub-MeV 
energy range. The HE emission also lasted much longer. 
The observations indicated a significan t devi ation from the 
so-called 'Band function' (|Band et all 1 19931 ') . namely, the 
presence of a hard power-law compon ent that dominates at 
high energies (e.g., JAbdo et al . 2009a) ). Ob servations of the 
short GRB 090510 (|Ackermann et alll20l61 ) were especially 
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fruitful: first, they proved the intrinsic similarity of the HE 
emission between short and long bursts; second, the detec- 
tion of a photon with energy over 30 GeV made it possible to 
evaluate the Lorentz factor of the jet F ~ 1000; finally, the 
simultaneous detection of signals in different energy ranges 
of GBM and LAT allowed one to put constraints on some 
quantum gravity theories which predict Lorent z invariance 
violation and energy-depend ent speed of light (|Abdo et al.l 
l2009bl : [GhTrlanda et alJl20ld ). 

In this letter we use the Fermi-LAT data to search for 
a high-energy emission related to GRBs. A similar question 
wa s addressed in a number of recent papers. In the papers 
by lAkerlof et all (|2010l . l201lT ) the matched filtering tech- 
nique was implemented and thr ee detections of HE emi ssion 
were claimed. In the paper by iBeniammT et al.1 i feOllT ) the 
HE fluence associated with several bright GBM events was 
constrained. 

Unlike the above studies, we looked for a HE emission 
that could extend over longer time spans than were exam- 
ined previously and which therefore could have been missed. 
For this purpose we searched the LAT database for HE pho- 
tons in two energy bands (E > 100 MeV and E > 1 GeV) 
and selected those which came from the directions of 581 
previously detected GRBs within the time window of 1500 s 
either before or after the burst. We then estimated the sta- 
tistical significance of the excess by comparing the observed 
number of photons with the expected background, taking 
into account the penalty for the total number of trials. The 
detection criterion was chosen in such a way that the entire 
search would give a single false detection with the probabil- 
ity of 0.05. We found 4 previously unreported GRBs showing 
the post-burst HE emission, and none showing the pre-burst 
emission. 



2 DATA 

The LAT detector has a sensitivity to photons with energy 
above 30 MeV with a wide field of view (FOV) of ~ 2.4 sr 
and the effective area of up to 9500 cm 2 . The detector an- 
gular resolution is a function of the incident photon energy 
and the 'event class' which is determined by the set of recon- 
struction cuts (Rando 2009). It also depends on the angle 
between the instrument axis and the arrival direction of a 
photon, but for the purpose of present analysis, we use the 
mission-average value. An extensive review of the instru - 
mental capabilities can be found in (|Atwood et al.ll2009T ). 
In this letter we use the LAT weekly all-sky data that are 
publicly available at Fermi mission websitqj. The analysis 
covers the time period of 127 weeks from August 04, 2008 
to January 07, 2011, corresponding to mission elapsed time 
(MET) from 239557417 s to 316111862 s. 

We use the 'diffuse' event class and impose an Earth 
relative zenith angle cut of 105°. We discard the photons 
from 'transient' and 'source' classes. The inclusion of these 
classes would decrease the signal-to-background ratio which 
becomes an impo rtant factor f o r our comparatively long 
time window (see lAkerlof et al.l (|201ll ) for detailed discus- 
sion of event classes with respect to GRBs). We do not re- 



quire the rocking angle cut of 52° in order to keep photons 
observed in the pointing mode, including repoint requests 
caused by the GRB trigger. 

For the spatial selection of photons we adopt the 95% 
containment angle ags(E,v) corresponding to th e point 
spread function (PSF) for the "diffuse' class photons (|Randol 
l2009l ; lAbdo et al.ll2009dl ; lBurnett et al.ll2009l ). This angle de- 
pends on both the photon energy and the conversion type 
v. The latter takes two discrete values: and 1 for front and 
back converted photons, respectively. 

In our analysis we used the times and coordinates of 
the GRBs detected by other instruments such a s GBIv^f] 
(iMeegan et al.ll2009h, SwiftR (iGehrels et alj|2004l), INTE - 
GRAlJ 4 ! (|Winkler et al.ll2003F. MAXI dMatsuoka et alj|2009h 
and Konus-Wind jAptekar et al.l 1 1993 ) ■ We have compiled 
two non-overlapping list of GRBs. The first list included 
605 GRBs detected by Fermi GBM and the second 279 
GRBs detected by other instruments only. Of these, 444 
and 137 GRBs, respectively, were in the Fermi-LAT FOV 
at least for some part of 1500 s after the burst. Note that 
although we used time stamps from the original observa- 
tions, localizations in many cases were provided by much 
more precise follow-up observations; these data were ob- 
tained through the Gamma-ray bursts Coordinates Network 
(GCN; |http://gcn.gsfc .iiasa.gov/ ). 



3 METHOD 

The key quantity in our analysis is the probability p that the 
observed HE emission from the direction of a given GRB is 
a fluctuation of the background. If this probability is smaller 
than the certain threshold, we claim the detection of the high 
energy emission from that burst. The significance threshold 
is obtained by requiring that the number of false detections 
does not exceed 0.05 in the entire set. Taking into account 
that the total number of bursts is 581 and counting two en- 
ergy ranges as independent, one obtains the following con- 
dition: 



p < 5 x 10~ 



(1) 



for either of the two energy regions. 

The probability p for a given burst and given energy 
threshold Eo is calculated as follows. Let t b , h and b b be the 
trigger time and galactic coordinates of a GRB. First, we de- 
termine the observed number of photons n above the energy 
-Eo by counting photons satisfying the following conditions: 



E>E , 



a(l, b, l b , b b ) < ^Jal 5 (E,v) + a 2 GKB , 
t b ^t^t b + 1500 s, 



(2) 



http: / / fermi.gsfc.nasa.gov/ssc / data/ access / 



where t, I, b, E and v stand for arrival time, coordinates, 
energy and conversion type of a photon, a(l,b,l b ,b b ) is the 
angular separation between photon and GRB, and ogrb is 
the GRB pointing error. The energy threshold Eo is either 
100 MeV or 1 GeV. These conditions select photons with 

2 http:/ /fcrmi. gsfc.nasa.gov/ssc/data/access/ 

3 http:/ /swift. gsfc.nasa.gov/docs/swift/ 

4 http:/ /www.isdc.unige.ch/integral/ 
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Table 1. List of Fermi- LAT GRBs showing extended high-energy emission. B is the expected background, n is the observed number of 
photons, p is the probability that the signal is fluctuation of the background, £ is the exposure, and t is total time duration within 1500 
s interval when the angle between boresight of the telescope and the GRB position was less than 65° . Previously unreported candidates 
are marked with the star. 
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Table 2. List of previously detected Fermi-LAT GRBs missed by our algorithm. 



energies larger than Eq that arrived within 1500 s after the 
burst from the region of interest (ROI). The latter is a cir- 
cle with the energy-dependent radius determined by the two 
contributions: the error of the photon arrival direction and 
the error of the GRB position. Usually the first contribution 
dominates. The error of the GRB position was taken to be 
equal to 1° in the case of GBM bursts and 0.5° in the case of 
bursts detected by Swift. Errors for all other bursts were de- 
termined individually from the GCN website. The observed 
pre-burst photons are selected by an obvious modification 
of the conditions ([2]). 

Next, we calculate the expected background B corre- 
sponding to the energy E > Eq. Since GRBs are exceptional 
events, for the background calculation we may use the pho- 
tons from the same spatial region for the entire duration 
of the mission. Thus, the background is given by the total 
number of photons observed from the ROI during the whole 
mission, multiplied by the ratio of the exposure correspond- 
ing to 1500 s after the burst to the total exposure of ROI. 
The calculation which takes into account the energy depen- 
dence of ROI is presented in the Appendix lAl 

Finally, having calculated the observed number of pho- 
tons n and the expected background B, the probability p 



for the GRB in question is calculated from the Poisson dis- 
tribution, 

p = P(B,n), 

where V(B, n) is the probability to observe n or more events 
at B expected. If this probability satisfies the condition |TJ 
for at least one of the two energy regions of interest, we 
have a detection and include the corresponding GRB in the 
detection list, Table [1] 



4 RESULTS AND CONCLUSIONS 

Applying the method of Sect. [3] to 581 GRB we have 
achieved 19 detections of the post-burst HE emission, 
of which 4 (namely, GRB 081009, GRB 090720B, GRB 
100728A and GRB 100911) were previously unreported. All 
detections correspond to GRBs present in the Fermi-GBM 
part of the GRB list. No pre-burst HE emission was found. 
Of the new detections, GRB 100728A demonstrated partic- 
ularly bright and long HE afterglow: 4 photons with energy 
> 1 GeV were observed vs. 0.065 expected from the back- 
ground, the chance probability being p = 7 ■ 10 -7 . Photon 
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arrival times relative to the burst trigger time are 711.3 s, 
713.8 s, 1161.0 s and 1342.5 s. GRB 081009, GRB 090720B 
and GRB 100911 were observed with the rocking angle less 
than 52°, while during the GRB 100728A the rocking an- 
gle exceeded 52° 1220 s after the burst. The low exposure 
for GRB100911 is because that burst went behind the Earth 
shortly after the trigger and there exists a possibility of con- 
tamination from the bright Earth limb. 

The total number of previously reported GRBs detected 
by Fermi-LAT is 20. Our algorithm failed in 5 cases; they 
are listed in Tabled Two of them (GRB 081215A and GRB 
100707A) resided far from the Fermi-LAT axis at approxi- 
mately 90° angular distance. Their original detections were 
made with the u se of the non-standard analys is technique 
(|McEnervl 120081 ; iPelassa fc Pesce-Rollinsj |201Ch . Our algo- 
rithm is based on the standard event reconstruction and 
does not treat events outside of LAT FOV. The remaining 
three GRBs are seen as excesses that do not satisfy eq. {]]) 
(in the worst case of GRB 100225 the chance probability 
is as large as 0.4). This discrepancy can be attributed to 
a wider time window used in our analysis and somewhat 
different energy ranges. 

In our procedure we have treated the background events 
as a Poissonian process. This approach would fail in case of 
a moving gamma ray source (the Sun or the Moon) cross- 
ing the region of interest just at the moment of burst. We 
have explicitly checked that no such events happened for the 
reported candidates. We have also assumed that the back- 
ground flux is stationary. This could lead to an erroneous 
detection if some gamma-ray sources in the region of interest 
flared at the moment of the GRB. This issue should be inves- 
tigated s eparately; one can use the LAT 1-year Point Source 
Catalog (|The Fermi-LAT Collaboration lioiiil ) for this pur- 
pose. No known sources appear within 95% containment ra- 
dius of all GeV photons attributed to GRB 100728A; the 
same is valid for GRB 081009. On the contrary, there are 
numerous sources in the neighborhoods of GRB 090720B 
and GRB 100911. The question of a possible influence of 
variability of thes e sources on the pre sent analysis will be 
studied elsewhere (|Rubtsov et aljfeoill ). Finally, to test for 
a possible influence of the magnetospheric flare^j we calcu- 
lated the gamma-ray flux during 1500 s after the burst in the 
ring between 15° and 20° from its location, and compared 
it to the expected background. No indication of magneto- 
spheric flares coincident with the reported new detections 
was found. 

When this paper was already submitted, analysis from 
the Fermi collaboration appeared confirming dete ction of 
HE emission from GRB 100728A (|Abdo et al.ll201ll ). 
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APPENDIX A: BACKGROUND CALCULATION 

The background expected from the region of interest dur- 
ing 1500 s after the burst is determined by the product of 
the differential flux f(E, I, b) from the ROI and the detector 
exposure in that direction. 

The number of photons coming within the time period 
[ti,ta] is: 

N = j dEdl dbf(E,l,b)£(E,l,b,ti,t 2 ), 

where £(E,l,b,ti,t2) is the exposure of the instrument in 
a given direction and the time period [ii,i2] at energy E, 
measured in units of cm 2 s. This exposure is estimated us- 
ing the standard Fermi-LAT tools gtltcube and gtexpeube 
(ScienceTools-v9rl8p6-fssc-20101108). 

The same number of photons can be written alterna- 
tively as 

N = J dEdl dbn(E,l,b,ti,t 2 ), 

where n(E, I, b, ti,t 2 ) is the angular spectral density of pho- 
tons detected over the period [ti , t 2 ] • Comparing these two 
equations one finds: 

n{E,l,b,ti,t 2 ) = f(E,l,b)£(E,l,b,ti,t 2 ). (Al) 

Therefore, knowing n[E,l,b,ti,t 2 ) and £(E,l,b,ti,t2) al- 
lows one to estimate the differential flux f(E, l,b). Under the 
assumption of a time-independent flux, the best estimate is 
obtained using the longest possible time period [i s tart , t cn d] . 

We represent n(E, I, b, istart, iend) as a discrete distribu- 
tion originated from all the photons observed for the whole 
duration of the mission [t start , Cendl ■ 

iV 

n(E,l,b,t stalt ,t cnd ) =J25(E-Ei) 8(l-h) 5{b-bi), (A2) 
»=i 

where h,bi,Ei are coordinates and energy of i— th photon. 
For the definition of the ROI we impose energy, spatial and 
time cut functions corresponding to Eq. 

a?«(E) = 6(E-E ), 

<T b (E, v, I, b) = 0(-a(l, b, l b , b b ) + sja% 5 (E,v) + a 2 GRB ) , 

a t {t) = 6(t - t b ) 6{-t + t b + 1500 s) , 

where 0(x) stands for Heaviside step function. 

Number of background photons is estimated as a sum of 
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front and back converted components each calculated with 
its own PSF and exposure: 

Be =J2[ dEdldbf{E,l,b) £(E, I, b, t b , t b + 1500 s) 



x af°(E) a b {E,v,l,b) 



(A3) 



where £ v stands for conversion-type dependent exposure 
(£ = £° + £ 1 ). We substitute / from EqUU and EqlAll 
to EalAH 



5(E - Ei) 8(1 - U) 5(6 - h) 



B E() = EE / dEdldh 



£(E,l,b,t start j ^cnd ) 



x £(E,l,b,t b ,t b + 1500 s) erf (E) a b (E, v, I, b) . 

Delta- functions will be integrated and E,l,b will be re- 
placed by Ei,h,bi. Sigma function will require to sum only 
the photons within ROI with energies larger than Eq. We 
finally get: 

N ^£ v (E % ,k,b l ,t b ,t b + 1500 s) 



q . ^ £ (-E'i j j 7 ^start j ^end ) 



where the sum formally goes through all LAT detected pho- 
tons, but effectively only ones in the spatial ROI do con- 
tribute. 
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